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Abstract 

Background: Powdery mildew (PM) is a major fungal disease of thousands of plant species, including many 
cultivated Rosaceae. PM pathogenesis is associated with up-regulation of MLO genes during early stages of 
infection, causing down-regulation of plant defense pathways. Specific members of the MLO gene family act as 
PM-susceptibility genes, as their loss-of-function mutations grant durable and broad-spectrum resistance. 

Results: We carried out a genome-wide characterization of the MLO gene family in apple, peach and strawberry, and 
we isolated apricot MLO homologs through a PCR-approach. Evolutionary relationships between MLO homologs were 
studied and syntenic blocks constructed. Homologs that are candidates for being PM susceptibility genes were inferred 
by phylogenetic relationships with functionally characterized MLO genes and, in apple, by monitoring their expression 
following inoculation with the PM causal pathogen Podosphaera leucotricha 

Conclusions: Genomic tools available for Rosaceae were exploited in order to characterize the MLO gene family. 
Candidate MLO susceptibility genes were identified. In follow-up studies it can be investigated whether silencing or a 
loss-of-function mutations in one or more of these candidate genes leads to PM resistance. 

Keywords: Rosaceae, MLO, Powdery mildew, Malus domestica 



Background 

Powdery mildew (PM) is a major fungal disease for thou- 
sands of plant species [1], including cultivated Rosaceae 
such as apple (Malus domestica), apricot {Prunus arme- 
niaca), peach {Prunus persica), and strawberry (Fragaria x 
ananassa). Powdery mildew occurs in all major growing 
regions of Rosaceous crops, leading to severe losses [2]. 
The major PM causal agents are Podosphaera leucotricha 
in apple [2], Sphaerotheca pannosa var. persicae in peach 
[3], Podosphaera tridactyla in apricot [4] and Podosphaera 
aphanis (syn, Sphaerotheca macularis f. sp. fragariae) in 
strawberry [5]. Powdery mildew shows similar symptoms 
in the four species: white spots appear on young green tis- 
sues, particularly leaves in the first days after opening. 
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whereas mature leaves show some resistance. Infected 
leaves crinkle, curl, and prematurely drop. Blossoms and 
fruits are not the primary targets of PM fungi, but infec- 
tions of these tissues are possible [2,3,5]. In peach, apricot 
and apple, PM spores overwinter in buds and then in 
spring, with the reprise of vegetative growth, the spores 
start a new infection [2,3]. 

Cultivars resistant to PM are fundamental in order to re- 
duce the use of pesticides in agricultural practices. The 
usual strategy in breeding focuses on dominant plant re- 
sistance genes (R-genes), however these genes often origin- 
ate from wild-relatives of the cultivated species, and thus 
interspecific crossability barriers could prevent their 
introgression [6]. Moreover, in case of a successful 
cross, several undesirable traits are incorporated with 
the R-gene, making extensive backcrossing necessary, 
which is time-consuming in woody species. Finally, the 
durability of R-genes is generally limited due to the ap- 
pearance of virulent strains of the pathogen, which can 
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overcome resistance in a few years [7] . Two examples are 
Venturia inaequalis race 6, which is able to overcome 
Rvi6 resistance to scab in apple [8], and P, leucotricha 
strains able to breakdown Pl-1 and Pl-2, two major PM 
R-genes of apple [9]. 

An alternative to the use of R-genes is based on plant 
susceptibility genes (S-genes), defined as genes whose 
loss-of-function results in recessively inherited resistance 
[10]. Barley mlo PM resistance, first characterized in 1942, 
is a remarkable example of immunity due to the absence 
of an S-gene, as it derives from a loss-of-function muta- 
tion in the MLO {Mildew Locus O) gene, encoding for a 
protein with seven transmembrane domains [11,12]. mlo 
resistance has long been considered as a unique form of 
immunity, characterized by durability, broad-spectrum 
effectiveness and recessive inheritance [13]. However, 
the characterization of the sources of resistance in other 
plant species, such as Arabidopsis [14], pea [15,16] and 
tomato [17], has revealed that resistance resulting from 
loss-of-function mutations in MLO functional orthologs 
is more common than previously thought. Therefore, it 
has been suggested that the inactivation of MLO sus- 
ceptibility genes could represent a valid strategy to 
introduce PM resistance across a broad range of culti- 
vated species [10]. 

Histological characterization of mlo resistance revealed 
that it is based on a pre-penetration defense system, asso- 
ciated with the formation of cell-wall appositions [14,18] 
and at least partially dependent on the actin cytoskeleton 
[19]. It has been suggested that functional MLO proteins 
negatively regulate vesicle-associated and actin-dependent 
defense pathways at the site of attempted PM penetration 
[20]. MLO proteins are therefore targeted by PM fungi as 
a strategy to induce pathogenesis. The early stages of PM 
infection are associated with an increase in transcript 
abundance of MLO susceptibility genes, showing a peak at 
six hours after inoculation. This has been shown to occur 
in tomato [17], barley [21], pepper [22] and grape [23,24]. 

MLO susceptibility genes are members of a gene family 
which shows tissue specific expression patterns and are in- 
volved in a variety of physiological processes, besides the 
response to PM fungi: one of the 15 MLO genes of Arabi- 
dopsis, AtML07, is involved in pollen tube reception by 
the embryo sac and its mutation results in reduced fertility 
[25]. Two other Arabidopsis genes, named AtML04 and 
AtMLOll, are involved in the control of root architecture, 
as mutants with null alleles of these two genes display 
asymmetrical root growth and exaggerated curvature [26]. 

Previous phylogenetic analysis of the MLO protein fam- 
ily identified six clades [23]. Currently, all MLO proteins 
functionally related to PM susceptibility in dicot species 
appear in a single clade, namely Clade V [14,17,23,24]. 
Similarly, Clade IV harbours all characterized PM suscep- 
tibility proteins from monocots [20,27]. 



MLO genes have been intensively studied in many 
monocots and dicots, but very little has been performed in 
Rosaceae. In this investigation, we characterized the MLO 
gene family in a number of Rosaceous species, with re- 
spect to their structural, genomic and evolutionary fea- 
tures. Moreover, we monitored the transcript abundances 
of apple MLO homologs following P, leucotricha inocula- 
tion in three apple cultivars. 

Results 

In silico and in vitro characterization of Rosaceae 
MLO homologs 

A database search for Rosaceae MLO homologs produced 
21 significant matches in peach, 23 in strawberry and 28 in 
apple. Of these, six (five from M. domestica and one from 
F. vesca) showed a very limited alignment region with other 
MLO genes, whereas eight (two from M. domestica, two 
from P, persica and four from F, vesca) were characterized 
by markedly different length with respect to MLO homo- 
logs reported in the genomes of Arabidopsis and grapevine 
[23,28], i.e. less than 350 amino acids (aa) or more than 700 
aa. Details on genomic localization amino acid number, pu- 
tative transmembrane domains and predicted exon/intron 
structure of the remaining homologs, together with infor- 
mation about the MLO homologs nomenclature chosen in 
this study is provided in Tables 1, 2 and 3. 

Peach and apricot are evolutionary very close to each 
other, and show a high degree of homology in DNA se- 
quence. Phylogenetic analysis (see next paragraph) indi- 
cated peach homologs PpMLOl, PpMLOS and PpML04 
as candidates for being required for PM susceptibility. 
Therefore, we used the sequences of these genes to design 
primers to identify full-length apricot MLO genes. This 
approach resulted in the amplification and the successive 
characterization of three MLO sequences, which were by 
analogy named PaMLOl, PaMLOS, and PaML04 (depos- 
ited in the NCBI database with the accession numbers 
KF177395, KF177396, and KF177397, respectively). 

Phylogenetic relations and inference of orthology 

We performed a phylogenetic study on the newly identi- 
fied Rosaceae MLO proteins. The dataset was completed 
with four homologs recently characterized in Rosa hybrida 
(rose) [30] (RhMLOl, RhML02, RhML03 and RhML04), 
the complete Arabidopsis thaliana AtMLO protein family 
[14], a series of MLO homologs which have been func- 
tionally associated with PM susceptibility, namely tomato 
{Solanum lycopersicum) SlMLOl [17], pea {Pisum sati- 
vum) PsMLOl [15,16], pepper {Capsicum annuum) 
CaML02 [27], lotus {Lotus japonicus) LjMLOl [15], barrel 
clover {Medicago truncatula) MtMLOl [15], barley {Hor- 
deum vulgare) HvMLO [11], rice {Oryza sativa) OsML02 
[31], wheat {Triticum aestivum) TaMLO_Bl and TaM- 
LO_Alb [31], and grapevine {Vitis vinifera) VvML014, 
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Table 1 Members of the MdMLO gene family as predicted in M. domestica cv. Golden Delicious genome sequence 



Gene Accession number^ Chr. Starting position (Mb) Clade Introns TM'' Amino acids Conserved aa*^ 



MdMLO 1 


MDPOOOOl 77099 


2 


1.02 


II 


11 


3 


487 


25 


MdML02 


MDP0000240125 


2 


11.10 


1 


11 


3 


571 


20 


MdMLOS 


MDPOOOOl 68575 


2 


11.11 


1 


13 


7 


670 


22 


MdML04 


MDP0000207002 


2 


8.79 


III 


16 


7 


634 


28 


MdMLOS 


MDPOOOOl 63089 


9 


15.26 


V 


14 


6 


579 


30 


MdML06 


MDPOOOOl 19433 


3 


33.95 


II 


0 


7 


504 


30 


MdMLOJ 


MDPOOOOl 23907 


n.d. 


n.d. 


V 


n.d. 


6 


561 


28 


MdMLOS 


MDP0000218520 


2 


11.11 


1 


9 


4 


390 


14 


MdML09 


MDP0000320797 


2 


27.20 


II 


10 


5 


454 


28 


MdMLO 10 


MDPOOOOl 96373 


3 


26.97 


1 


13 


6 


539 


28 


MdMLO n 


MDP0000239643 


4 


9.84 


V 


12 


8 


575 


28 


MdMLO 12 


MDPOOOOl 331 62 


6 


0.81 


III 


13 


5 


516 


28 


MdMLOlS 


MDPOOOOl 42608 


7 


7.48 


II 


12 


6 


351 


18 


MdMLO 14 


MDPOOOOl 91 469 


8 


29.25 


II 


10 


5 


395 


23 


MdMLOlS 


MDPOOOOl 41 595 


9 


7.54 


III 


15 


6 


647 


24 


MdMLOie 


MDPOOOOl 91 848 


9 


21.12 


VI 


14 


6 


606 


29 


MdMLO 17 


MDPOOOOl 45097 


11 


27.97 


1 


13 


7 


523 


28 


MdMLOlS 


MDP0000928368 


10 


27.97 


VII 


12 


7 


502 


30 


MdMLO 19 


MDPOOOOl 6871 4 


12 


16.23 


V 


13 


7 


590 


30 


MdMLO20 


MDPOOOOl 34649 


13 


11.61 


VIII 


13 


5 


589 


27 


MdML021 


MDPOOOOl 33760 


15 


24.99 


VI 


15 


6 


560 


28 



^Available at http://www.rosaceae.org/gb/gbrowse/malus_x_domestica/. 

"^Number of transmembrane domains in the predicted protein, as determined by InterPro prediction software (http://www.ebi.ac.uk/interpro/). 
'^number of conserved amino acids out of the 30 identified by Elliot et al. [29]. 



the only dicot MLO homolog known to belong to clade IV 
[23]. Clustering analysis using the UPGMA algorithm re- 
sulted in a total of eight distinct clades and no divergent 
lineage (Figure 1). Clade numbers from I to VI were 
assigned based on the position of Arabidopsis AtMLO ho- 
mologs and barley HvMLO, according to the previous 
study of Feechan et al [23]. The two additional clades 
(named VII and VIII) were found to include Rosaceae 
MLO homologs only, both having one homolog from P, 
persica, one from F, vesca and one from M, domestica. Fur- 
ther clustering analysis with a Neighbour- Joining algorithm 
resulted in merging clade VII and VIII (not shown). 

Four apple MLO homologs (MdMLOS, MdML07, 
MdMLOll and MdML019) and three MLO homologs 
from peach (PpMLOl, PpML03 and PpML04), apricot 
(PaMLOl, PaML03 and PaML04) and woodland straw- 
berrry (FvMLOl, FvML04 and FvML012) were found 
to cluster together in the phylogenetic clade V, contain- 
ing all the dicot MLO proteins experimentally shown to 
be required for PM susceptibility (e.g. [16,23]. One 
homolog from strawberry (FvMLOl 7) and one from 
peach (PpMLOl 2) were found to group, together with 
grapevine VvML014, in clade IV, which contains all 



monocot MLO proteins acting as PM susceptibility fac- 
tors (Figure 1). 

We used the GBrowse-Syn tool to detect syntenic blocks 
encompassing P, persica, F. vesca and M, domestica MLO 
genes. As syntenic blocks derive from the evolution of the 
same chromosomal region after speciation, orthology be- 
tween MLO genes could be inferred. In total, twelve ortho- 
logous relationships were predicted between P, persica and 
F, vesca, nine between P, persica and M, domestica and 
eight between F, vesca and M. domestica (Table 4, Figure 2 
and Additional file 1). 

Transcription of putative apple MLO genes in response to 
Podosphaera leucotricha inoculation 

To identify MLO genes that respond to the PM fungus P. 
leucotricha, we measured the transcript abundance of 19 
out of 21 apple MLO genes in leaves 4, 6, 8 and 24 hours 
after artificial inoculation with the pathogen, and compared 
these data with the ones of non-inoculated leaves. Three 
cultivars. Golden Delicious, Braeburn and Gala, were ana- 
lysed in order to investigate whether up-regulation was 
comparable among them and results could therefore be 
generalized for all apple cultivars. Three genes, namely 
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Table 2 Members of the PpMLO gene family as predicted In Prunus persica genome sequence 


Gene 


Accession number^ 


Chr. 


Starting position (Mb) Clade 


Introns 


tm'' 


Amino acids 


Conserved aa*^ 


PpMLO 1 


ppa003207m 


6 


6.82 V 


14 


7 


593 




30 


PpML02 


ppa003435m 


7 


18.38 III 


14 


8 


574 




30 


PpMLOS 


ppa003437m 


6 


21.99 V 


13 


7 


574 




30 


PpML04 


ppa003466m 


2 


21.03 V 


14 


7 


572 




30 


PpMLOS 


ppa003706m 


4 


10.92 1 


14 


8 


555 




30 


PpMLOe 


ppa004012m 


7 


22.64 II 


14 


6 


535 




29 


PpMLOJ 


ppa004508m 


8 


21.17 II 


0 


7 


506 




29 


PpMLOS 


ppa004621nn 


6 


22.01 VI 


14 


6 


499 




29 


PpML09 


ppa004687nn 


4 


2.59 VII 


11 


7 


496 




29 


PpMLO 10 


ppa004866m 


2 


13.73 II 


11 


7 


488 




29 


PpMLOll 


ppa020172m 


1 


43.04 1 


14 


4 


561 




30 


PpMLOU 


ppa020311m 


5 


0.82 IV 


13 


7 


566 




30 


PpMLOlS 


ppa02 1048m 


4 


15.57 VIII 


12 


5 


510 




24 


PpMLOH 


ppa022847m 


6 


6.80 VI 


14 


6 


550 




29 


PpMLOlS 


ppa024476m 


7 


1 7.63 1 


14 


8 


539 




26 


PpML016 


ppa024488m 


5 


0.76 III 


14 


6 


504 




30 


PpMLOU 


ppa024581m 


6 


8.95 II 


13 


6 


463 




27 


PpMLOlS 


ppa026565nn 


6 


22.00 VI 


13 


6 


416 




25 


PpMLO 19 


ppb024523m 


1 


42.04 II 


13 


5 


446 




23 


^Available at http://www.rosaceae.org/gb/gbrowse/prunus_persica/. 

"^Number of transmembrane domains in the predicted protein, as determined by InterPro prediction software (http://www.ebi.ac.uk/interpro/). 
^number of conserved amino acids out of the 30 identified by Elliot et al. [29]. 






Table 3 Members of the FvMLO gene family as predicted In Fragaria vesca genome sequence 






Gene 


Accession number^ 


Chr. 


Starting position (Mb) Clade 


Introns 


tm'' 


Amino acids 


Conserved aa*^ 


FvMLO 1 


mrna02774.1-vl.0-hybrid 


n.d. 


n.d. V 


14 


7 


632 


28 




FvML02 


mrna03210.1-vl.0-hybrid 


3 


14.46 II 


11 


5 


528 


20 




FvMLOS 


nnrna09651.1-vl.0-hybrid 


6 


35.88 III 


14 


6 


542 


28 




FvML04 


mrna09653.1-vl.0-hybrid 


6 


35.90 V 


14 


7 


573 


30 




FvMLOS 


mrnal0166.1-vl.0-hybrid 


1 


1.34 II 


14 


3 


688 


26 




FvML06 


mrnal0346.1-vl.0-hybrid 


3 


12.52 II 


7 


2 


385 


15 




FvMLOJ 


mrnal0363.1-vl.0-hybrid 


3 


12.49 II 


9 


2 


442 


21 




FvMLOS 


mrnal 0558.1 -vl.O-hybrid 


2 


19.08 II 


n.d. 


6 


514 


28 




FvML09 


mrnall028.1-vl.0-hybrid 


n.d. 


n.a. 1 


10 


4 


434 


18 




FvMLO 10 


nnrnal3023.1-vl.0-hybrid 


1 


7.96 III 


13 


6 


557 


27 




FvMLO 11 


mrnal 4592.1 -vl.O-hybrid 


1 


8.77 1 


13 


7 


548 


28 




FvMLO 12 


mrna231 98.1 -vl.O-hybrid 


7 


15.89 V 


14 


7 


507 


29 




FvMLOlS 


mrna26428.1 -vl.O-hybrid 


7 


17,79 VIII 


11 


5 


558 


20 




FvMLOH 


mrna28541.1 -vl.O-hybrid 


n.d. 


n.a. Ill 


11 


4 


481 


26 




FvMLO IS 


mrna29770.1 -vl.O-hybrid 


3 


7.36 VII 


13 


7 


538 


28 




FvMLOie 


mrna31264.1-vl.0-hybrid 


3 


30.51 1 


16 


8 


579 


28 




FvMLO 17 


mrna31 498.1 -vl.O-hybrid 


5 


20.23 IV 


11 


5 


531 


27 




FvMLO IS 


mrna29285.1 -vl.O-hybrid 


5 


19.12 V 


6 


4 


357 


18 





^Available at http://www.rosaceae.org/gb/gbrowse/fragaria_vesca_v1.0-lg/ (hybrid). 

"^Number of transmembrane domains in the predicted protein, as determined by InterPro prediction software (http://www.ebi.ac.uk/interpro/). 
^number of conserved amino acids out of the 30 identified by Elliot et al. [29]. 
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Figure 1 Phylogenetic tree of Rosaceae MLO. Phylogenetic relationships of predicted Rosaceae MLO amino acid sequences to MLO proteins 
of other plant species. The dataset includes Rosaceae MLO sequences from Roso hybrido (RhMLO), Molus domestica (MdMLO), Prunus persica 
(PpMLO), Prunus ormenioco (PaMLO) and Frogaria vesco (FvMLO). The other proteins included are Solonum lycopersicum SMLOI , Arabidopsis thaliono 
AtMLO, Capsicum annuum CaML02, Pisum sativum Psf\AlO] , Medicago truncatula MtMLOl, Lotus japonicus LjMLOl, Vitis vinifera VvML014, Hordeum 
vulgare HvMLO, Triticum aestivum TaMLO_Bl, TaMLO_Alb and Oryza sativa OsML02. Proteins which have been functionally characterized as 
susceptibility genes are highlighted in bold. Numbers at each node represent bootstrap support values (out of 100 replicates). 



MdMLOlh MdMLOlS and MdML019, were found to 
be significantiy up-regulated after inoculation with the 
pathogen (Figure 3 and Additional file 2). Up-regulation of 
these genes was about 2-fold compared to non-inoculated 
plants, with peaks of 4-fold up-regulation at very early time 
points CBraeburn - MdMLOll - 6 hpi; 'Gala - MdMLOlS - 
4 hpi; 'Golden Delicious'- MdML019 - 6hpi), MdMLOll 
and MdMLOlS were up-regulated in all cultivars, while 
MdMLO 19 was only up-regulated in 'Braeburn' and 
'Golden Delicious'. 

Two of the genes, MdMLOll and MdMLO 19 belong 
to Clade V, while MdMLOlS belongs to the newly iden- 
tified Clade VII (Figure 1), 



Discussion 

Genomic organization and phylogenetic relations 
between Rosaceae MLO homologs 

We report here the identification, through an in silico ap- 
proach, of 19 MLO homologs in the genome of peach and 
18 in the genome of strawberry. This is consistent with the 
results of previous genome-wide studies carried out on di- 
cotyledonous species, indicating the presence of 15 MLO ho- 
mologs in Arabidopsis, 17 in grapevine and 16 in tomato 
[9,13]; Appiano et at, unpublished results; [24]. Conversely, 
the number of MLO homologs detected in apple (21) is lower 
than expected, considering that a relatively recent genome- 
wide duplication event has occurred in the Pyreae tribe [32]. 
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Table 4 Relations of orthology Inferred between P. persica, 
F, vesca and M. domestica MLO homologs 



P. persica genes F. vesca orthologs M. domestica orthologs 



PpML02 


FvMLOlO 


MdMLOlS 


PpMLOS 


FvML04 


MdML019 


PpML04 


FvMLOU 


- 


PpMLOS 


FvMLOW 


MdMLOlO, MdMLOU 


PpMLOe 


FvMLOS 


MdMLOl 


PpMLO? 


FvML08 




PpMLOS 


FvMLOS 




PpML09 


FvMLOlS 


MdMLOlS 


PpMLO 10 


FvML02 


MdML09 


PpMLOH 


FvMLOM 


MdML021 


PpMLOlS 


FvMLOU 




PpMLO 16 




MdMLOU 


PpMLO 18 


FvMLOS 





Relations of ortliology between PpMLOl, PpMLOS, PpML04 and apricot 
PaMLOl, PaMLOS, PaML04 were clearly suggested by the high percentage of 
sequence identity between these homolog genes, which was 97,3%, 98,8% 
and 96,7%, respectively. 



Most PpMLO, FvMLO and MdMLO homologs appeared 
to be widely distributed within the respective genomes 
(Tables 1, 2 and 3), indicating segmental duplication as the 
prevailing evolutionary mechanism for the Rosaceae MLO 
gene family. However, we also found cases of clusters of 
adjacent homologs {PpMLOS, PpMLOS and PpMLOlS, 
PpML012 and PpMLOlG, PpMLOl and PpML014, 
FvMLOS and FvML04, FvML06 and FvML07, MdML02, 
MdMLOS and MdMLOS), which are likely the result of 
tandem duplication events. 

Inference of phylogenetic relationships between MLO 
proteins revealed the presence of apple, strawberry, peach 
and apricot homologs in the clade V, containing all dicot 
MLO homologs shown so far to be involved in PM sus- 
ceptibility, thus making them candidates to act as suscep- 
tibility factors. Although the simple clustering in clade V 
is not enough to recognize a gene as a susceptibility factor, 
it does provide the first evidence for functionality and al- 
lows for the reduction in the number of candidates for fur- 
ther functional analysis. Clade IV, that contains functional 
MLO susceptibility homologs from monocots, was found 
to include one homolog from F. vesca (FvMLO 17) and 
one from P, persica (PpMLO 12). In accordance with this 
finding, a MLO homolog from the dicot species V. vinifera 
also clusters in clade IV [23,24]. Figure 1). Interestingly, 
phylogenetic analyses carried out in this study also re- 
vealed the presence of one or two additional clades, 
depending on the type of phylogenetic reconstruction 
(UPGMA or Neighbour- Joining), which were not reported 
to occur in earlier investigations. Moreover, they appear to 
be characteristic of Rosaceae, since they contain only 



homologs from this family. Clearly, the exclusivity for 
Rosaceae of these clade(s) needs to be confirmed by fur- 
ther studies containing a larger dataset of MLO pro- 
teins. Additional studies could be also addressed to the 
functional characterization of Rosaceae MLO homologs 
grouped in clade VII. Indeed, this appears to be basal to 
both clade IV and clade V (Figure 1), and thus might have 
contained ancestral proteins which later on evolved into 
PM susceptibility factors. 

Synteny between apple, peach and woodland strawberry 
MLO genes 

We took advantage of recent developments in Rosaceae 
genomics in order to detect synteny between P, persica, 
F, vesca and M, domestica chromosomal regions con- 
taining MLO homologs. This permitted the inference of 
ortholgous relationships between MLO genes in these 
species. Notably, all predicted MLO orthologs from dif- 
ferent Rosaceae species, fell in the same phylogenetic 
clade (Tables 1, 2 and 3; Figure 1 and Additional file 1). 
This is to be expected, since orthologs generally share 
the same function and thus are characterized by a high 
level of sequence conservation. It is noteworthy that the 
chromosomal localization of predicted MLO orthologs 
between P, persica, M, domestica and F, vesca is in ac- 
cordance with the results of the synteny study performed 
after the release of the three genomes [33,34] . In particu- 
lar, genes situated on peach scaffold 2, 7 and 8 were pre- 
dicted to have orthologs on strawberry chromosome 7, 1 
and 2, respectively, whereas genes on peach scaffold 4 
were predicted to have orthologs on strawberry chromo- 
somes 2 or 3 (Table 4). FvMLOS was predicted to be 
orthologous to two peach MLO genes, PpMLOS and 
PpMLO IS, which are localised in close proximity to each 
other on peach scaffold 6 and grouped together in clade 
VI. In this case, we hypothesize a relation of co-orthology 
due to the occurrence of a recent tandem duplication 
event in the peach genome. Similarly, PpMLOS and 
FvMLO 16 were predicted to be orthologs of two apple 
MLO genes, MdMLO 10 and MdMLO 17, located on 
chromosomes 3 and 11. This is consistent with indica- 
tions of duplications of large segments of these two 
chromosomes during the evolution of the apple gen- 
ome [32]. 

Transcription of apple putative MLO genes In response to 
R leucotricha Inoculation 

In barley, pea and tomato, only one of the clade V MLO 
homologs seems to be involved in powdery mildew sus- 
ceptibility, whereas in A. thaliana three MLO genes in 
Clade V are required to be inactivated in order to achieve 
a fully resistant phenotype [16,27]. This implies that, 
within Clade V MLO genes, a further selection might be 
required to identify PM susceptibility genes. Accumulating 
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evidence indicates that MLO susceptibility genes are 
up-regulated upon challenge with powdery mildew 
fungi [17]. Therefore, we analysed the expression level 
of apple MLO genes identified in this study in response 
to the interaction with P. leucotricha. Three pathogen- 
dependent gene up-regulations were detected. Two up- 
regulated MLO homologs, MdMLOll and MdML019, 
encode for proteins falling in clade V, thus making 
them likely candidates to act as PM susceptibiUty genes 
in apple. MdMLOll and MdML019 are located on 
chromosomes 4 and 12 respectively, and are therefore 
both generated from a duplication event in the 9- 



chromosome ancestor of apple [32]. A third pathogen- 
dependent up-regulated gene, MdMLOlS, was found, 
which encodes a protein grouping in the newly identi- 
fied Clade VII (Figure 1). The presence of a PM upregulated 
gene outside clade V is consistent with transcriptome 
analyses recently performed in tomato (Appiano et al, 
unpublished results). Apple clade V also contains two 
genes, MdMLOS and MdML07, which show no sig- 
nificant changes in expression following inoculation. 
Accordingly, the lack of up-regulation of some clade V 
MLO genes has been observed in grapevine and tomato 
[23,24]; Appiano et al, unpubUshed results). The 



Pessina et al. BMC Genomics 2014, 15:618 
http://www.bionnedcentral.conn/1471 -21 64/1 5/61 8 



Page 8 of 1 2 



MdMLOll 



MdMLOlS 



MdML019 




2,5 
2 

1.5 
1 

0.5 
0 



2 
1.5 
1 

0.5 
0 

7 
6 
5 
4 
3 
2 
1 
0 



3.5 
3 

2.5 
2 



mm 'mm 



4h 6h 8h 24h 



4h 6h 8h 24h 



* * 










T 






1 


1 



4h 6h 8h 24h 



I — — r 

4h 6h 8h 24h 



4h 6h 8h 24h 




6h 8h 24h 



■ Control 



Significantly different values (p < 0.05): ^ up-regulation 

□ Inoculated 

Figure 3 Transcriptional variation of three apple MLO genes following inoculation with P. leucotricha. Transcript abundances of three 
MLO genes in leaves of the apple cultivars 'Braeburn', 'Golden Delicious' and 'Gala', following powdery mildew (PM) inoculation. Here we show 
only MLO genes that were more than one time significantly up or down regulated following PM inoculation at one of the four time points 
examined (4, 6, 8 and 24 hpi). The set of results of all investigated genes is shown in Additional file 1 : Figure SI . Each bar shows the average of four to 
eight biological replicates. The Ct values have been normalized for three reference genes: actin, ubiquitin and elongation factor 1. Statistical significance 
was determined with a t-test for each individual pair of inoculated and non-inoculated samples at each time point. The error bars show standard errors 
of the means. Significant differences between inoculated samples and control samples are indicated with an asterisk (P< 0.05). 



possible role of these genes as susceptibility factors has 
not yet been highlighted. 

PpMLOS, PaMLOS and FvML04 are likely to represent 
true orthologs of MdML019 (Table 4). Since orthologs 
often maintain the same function during evolution, we 
conjecture that the expression of these genes might 
also be responsive to PM fungi attacking correspond- 
ing species. Moreover, FvMLOlS and PpML09 are 
likely orthologs of MdMLOlS, so they should also 
be considered as putative transcriptionally responsive 
genes to PM fungi attack. Further studies aimed at the 
functional characterization of these genes (e.g. through the 
application of reverse genetic approaches of targeted mu- 
tagenesis or gene silencing), in apple but also in peach and 
strawberry, might lead to the identification of resistant 
phenotypes, which could be used for the development of 
PM resistant cultivars. Particularly, studies on MdMLOlS 



could lead to the characterization of a possible role for 
clade VII in the interaction with PM fungi. 

Conclusions 

Our work led to the identification of 19 MLO homo- 
logs in peach, 17 in strawberry and 21 in apple. Three, 
three and four homologs, respectively, belong to clade 
V and therefore are candidates for being S-genes. Due 
to the high similarity between peach and apricot, we 
were able to amplify and characterize three Clade V 
apricot MLO genes. 

The phylogenetic analysis revealed two new Rosaceae 
specific clades for the MLO family, although this needs 
to be confirmed by the use of a larger MLO proteins 
dataset. 

Through inoculation of apple with P. leucotrica, we 
identified three up-regulated genes, i.e. MdMLOll , 
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MdMLOlS and MdML019. MdMLOll and MdML019, 
that belong to Clade V, are positioned in duplicated re- 
gions and have high sequence identity, therefore they are 
likely to be recent paralogs. MdMLOlS belongs to the 
newly identified Clade VIL 

Methods 

In silico identification and comparison of MLO predicted 
proteins in peach, woodland strawberry and apple 

Predicted peptides from the peach genome (v. 1.0) and the 
strawberry genome (v. 1.0) gene model databases, available 
at the website of the Genomic Database for Rosaceae 
(www.rosaceae.org) [35], were queried for the presence of 
MLO homologs protein sequences. First, a BLAST search, 
using the tomato SlMLOl amino acid sequence as query 
was carried out. A further search was performed with the 
HMMER programme, which uses a method for homolog 
searches based on the profile hidden Markov probabilistic 
model [36]. The sequences obtained with the previously 
mentioned BLAST search, were used together with other 
known MLO sequences from dicot and monocot species, 
namely: four RhMLOs from Rosa hybrida, 15 AtMLOs 
from Arabidopsis thaliana, SlMLOl from Solanum lyco- 
persicum, CaML02 from Capsicum annuum, PsMLOl 
from Pisum sativum, MtMLOl from Medicago trunca- 
tula, LjMLOl from Lotus japonicus, VvML014 from V, 
Vinifera, HvMLO from Hordeum vulgare, TaML01_Alb 
and TaMLO_Bl from Triticum aestivum and OsML02 
from Oryza sativa, MLO protein sequences from apple 
{Malus domestica cv. Golden Delicious) were identified by 
searching for the MLO domain profile (IPR004326) in the 
apple genome available at FEM-IASMA computational 
biology web resources (http://genomics.research.iasma.it). 
The resulting list was integrated with a BLAST search, car- 
ried out with the amino acid sequences previously listed for 
the HMMER search in peach and strawberry. 

Chromosomal localization and predicted introns/exons 
structure of each MLO gene of apple, peach and strawberry 
was deduced based on the available genomic information at 
the GDR database. The presence and number of membrane 
spanning helices was predicted using the online software 
InterPro (http://www.ebi.ac.uk/interpro). Alignments for 
conserved amino-acids analysis were performed with the 
CLC Sequence Viewer v. 6.9 software (http://clcbio.com). 

A total of 90 MLO protein sequences, including three 
apricot MLO sequences isolated in vitro (see next para- 
graph), were used for Clustal alignment (http://www.ebi.ac. 
uk/Tools/msa/clustalw2/). UPGMA-based and Neighbour- 
joining-based phylogenetic trees were obtained with the 
CLC sequence viewer software. The UPGMA clustering al- 
gorithm was further used as input for the Dendroscope 
software, suitable for the visualization of large phylogenetic 
trees [37]. 



Relationships of orthology between MLO candidate 
genes from peach, woodland strawberry and apple 
were inferred by running the GBrowse-Syn tool avail- 
able at GDR (http://www.rosaceae.org/gb/gbrowse_syn/ 
peach_apple_strawberry) [35,38]. This displays syntenic 
regions among the three available genomes of Rosaceae, 
as detected by the Mercator programme [35,39]. For 50 
Kb chromosomal stretches flanking each P, persica 
PpMLO homolog, syntenic regions from F. vesca and M. 
domestica were searched. Orthology was called upon 
the identification of F. vesca or M. domestica MLO ho- 
mologs within syntenic blocks. 

In vitro isolation of apricot MLO homologs 

RNA from apricot leaves (cultivar Orange Red) was ex- 
tracted by using the SV Total RNA Isolation System Kit 
(Promega), and corresponding cDNA was synthesized 
by using the QuantiTect Reverse Transcription Kit (Qia- 
gen) with oligo(dT) primers. Sequences of the peach MLO 
homologs PpMLOl, PpMLOS and PpML04, are phylo- 
genetically close to MLO homologs functionally associated 
to PM susceptibility, and were therefore used to design 
the primer pairs 5'-ATGGCAGCCGCAACCTCAGG 
AAGA-3 75 '-TTATATACTTTGCCTATTGTCAAAC-3 

5 '-ATGGCAGGGGGAAAAGAAGGACG-3 75 -TCAAC 
TCCTTTCTGATTTCTCAA-3' and 5'-ATGGCCGA 
ACTAAGTAAAG A-3 75 ' TC AACTTCTTG ATTTTCC 
TTTGC-3', respectively. These were employed to 
amplify full-length cDNA sequences of apricot putative 
orthologs, by using the AccuPrime Taq polymerase 
(Invitrogen). Amplicons were purified by using the 
NucleoSpin Extract II kit (Macherey-Nagel) and ligated 
(molar ratio 1:1) into the pGEM-T Easy vector (Pro- 
mega). Recombinant plasmids were cloned in E. coli 
DHlOp chemically competent cells and recovered by 
using the Qiaprep spin miniprep kit (Qiagen). Sequen- 
cing reactions were performed twice, by using univer- 
sal T7 and SP6 primers (Eurofins MWG Operon). 

Glasshouse test with apple cultivars 

A total of 192 apple plants from three cultivars (Braeburn, 
Golden Delicious and Gala) were used to measure tran- 
script abundance of MLO genes. Budwoods from these cul- 
tivars were grafted on M9 rootstocks in January 2012. The 
grafts were kept at -1°C for 2 months, and potted at the be- 
ginning of March in greenhouse. The plants grew for 

6 weeks in the greenhouse at 20°C during the day, 17°C 
during the night, relative humidity of 70% and natural day/ 
night cycle. 

P. leucothrica was collected from apple trees in an un- 
sprayed test orchard and used to infect greenhouse grown 
apple seedlings from 'Gala Galaxy' seeds. Four weeks after 
inoculation, conidia were used for the inoculation experi- 
ment, or transferred to new seedlings, to keep them viable. 
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We inoculated by touching the plants with heavily infected 
apple seedlings. Control plants were not inoculated and 
kept separated in the same greenhouse of the inoculated 
plants. Inoculated and control plants were grown in the 
greenhouse at the growing conditions previously men- 
tioned. The leaf samples were collected 4, 6, 8 and 
24 hours post-inoculation (hpi). 

Eight experimental repeats were performed and each 
sample contained three or four young leaves collected 
from each single plant. Every plant was used for sampling 
only once, to avoid any possible effect of wounding on the 
expression ofMLO genes. The smallest statistical unit was 
a plant. The leaves were flash-frozen and ground in liquid 
nitrogen, and stored at -80°C until RNA extraction. 

qPCR analysis of transcript levels 

RNA extraction was carried out with the MagMAX-96 
Total RNA isolation kit (Applied Biosystem) that in- 
cludes DNAse treatment. The kit yielded between 50 
and 200 ng/ul, of good quality RNA per sample. 

Primers for gene expression analysis were designed with 
NCBI Primer Designing Tool (http://www.ncbi.nlm.nih. 
gov/ tools/primer-blast/). Four serial dilutions of cDNA (1/ 
5 - 1/25 - 1/125 - 1/625) were used to calculate the 



efficiency of each primer pair with iCycler software (Biorad). 
In case of efficiency lower than 1.80 or greater than 2.20, the 
primer pair was discarded and a new one tested, with the ex- 
ception of MdML09, for which was not possible to design a 
primer pair with better efficiency. It was only possible to ana- 
lyse 19 MLO genes because for MdML012 and MdML016 
was not possible to design specific and efficient primer pairs, 
despite numerous attempts. Presence of a specific final dis- 
sociation curve was determined after each qPCR amplifica- 
tion reaction with progressive increment of temperature 
fi:om 65°C to 95°C (0.5°C each step, 5 sec) and the size of the 
product was confirmed by agarose gel electrophoresis. 

Quantitative Real Time-PCR (qPCR) was performed 
with SYBR greenER mix (Invitrogen) in a 15-(iL reac- 
tion volume, using a Bio-Rad iCycler iQ detection sys- 
tem, run by the Bio-Rad iCycler iQ multicolor 3.1 
software. The software applies comparative quantifica- 
tion with an adaptive baseline. Samples were run in two 
technical replicates with the following thermal cycling 
parameters: 95°C 3 min - 95°C 15 sec, 60°C 1 min (re- 
peated 40 times) - 95°C 10 sec. 

Reference genes |3-actin (NCBI accession number 
DT002474; Plaza accession number MD00G171330 - 
http://bioinformatics.psb.ugent.be/plaza/), ubiquitin (Plaza 



Table 5 Gene-specific 
domestica cv. Golden 



primers and amplicon sizes in qRT-PCR detection of 19 MdMLO-Wke genes based on Malus 
Delicious genome sequence 



Gene 


Forward primer (5' - 3') 


Reverse primer (5' - 3') 


Size (bp) 


Efficiency 


MdMLOl 


GTGGGCTCGGTCGGCCAAAA 


CCAGCACCAGCACCAGAACCA 


81 


2.06 


MdML02 


CGTOGATCAACCACTGCGCCT 


TGAGCTGCAGCCAGTGGGATCT 


87 


1.83 


MdMLOS 


CCACTGCGCCTCTCTGAAGCA 


CCACCAAAACGGCTCTCCAGGT 


93 


2.12 


MdML04 


TGTOCAGACACTATGCTGCCATGT 


GGCAGCAGCTAAAGATCTGCGT 


109 


1.87 


MdMLOS 


TCGTCAGGCTCTCATOGGGGT 


GTGCTGCTGCCACTCCCTC 


132 


1.80 


MdML06 


TCGCGGAGGAGGGGTCG^ 


TOGAGCGACAGCAACGGCA 


72 


2.15 


MdML07 


TGGAGCAAGTCACCAGTCTCCAT 


CGOTCCTGGTGCCAAATGTGC 


127 


2.12 


MdMLOS 


GTCAAGCTAATCTOCCACGCGCT 


GGCTGGAAGGAAGGACAGCCA 


85 


1.95 


MdML09 


GCTGCAACACGTAATCACCC 


AGAACGCCAmCGAAAGCA 


173 


2.30 


MdMLOl 0 


GCGATCGTOGCOTGACTCC 


TOCGCGCTCGACAAGCAGA 


86 


1.92 


MdMLOll 


CCGTOCATCACCAAGACGA 


ATOCTCTCCGAGTOCGCC 


102 


1.90 


MdMLOl 3 


ACATOTCCCCAGGOTG^ 


GCCCAACCAATAAGTCCCGA 


151 


2.00 


MdMLOl 4 


TGCAOTGTCAGCCAGATGGG 


GCATCTCCCACCCACGAACCG 


81 


2.15 


MdMLOl 5 


GCGCCmCTCTCTGCTGGGT 


CGCGTGCGAGGTGGTCTOT 


90 


2.01 


MdMLOl/ 


TOCCACTGTATGCOTGGT 


TGOTGOTCTGTGCGAATG 


163 


2.15 


MdMLOl 8 


AAGGAAGGCTCTCATOAGGCTCT 


TGCAATOGC^GACCAACGGT 


100 


2.22 


MdMLOl 9 


CAGAGTGGCGACTGCAC™ 


GGGACATGGAGTGCAAAGGA 


110 


1.97 


MdMLOIO 


AAAAAGCTCCACCAACCCCA 


mCTCTCCCATGACGCTCG 


165 


2.11 


MdML021 


COTGTOGAGGCCGTAGAG 


ACCAAGTGCmGGTGGm 


176 


1.95 


P-actin 


CTATGTOCCTGGTATOCAGACC 


GCCACAACOTGA^CATGC 


82 


1.90 


Ubiquitin 


CATCCCCCCAGACCAGCAGA 


ACCACGGAGACGAAGCACCAA 


349 


1.91 


Elongation Factor 1 


TACTGGAACATCACAGGCTGAC 


TGGACCTCTCAATCATGTOTC 


308 


2.07 
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accession number MD05G001920) and elongation factor 
1 (Plaza accession number MD09G0 14760) were used as 
reference genes (Table 5). All these three genes were used 
in previous works [40-42]. For additional control, we 
assessed the stability of our genes with the software geN- 
orm (medgen.ugent.be/~jvdesomp/genorm/). An M-value 
lower than 1.5 is generally considered as stable enough 
[43-45] and all three reference genes in all three cultivars 
considered are within this threshold. We saw differences 
in stability between cultivars: 'Golden Delicious' was the 
most stable cultivar (actin: 0.824 - ubiquitin: 0,852 - 
elongation factor 1: 0,926), whereas 'Braeburn was the less 
stable (actin: 1.246 - ubiquitin: 1,293 - elongation factor 
1: 1,369) and 'Gala showed intermediate stability (actin: 
1.039 - ubiquitin: 1,152 - elongation factor 1: 1,078). 

Each of the biological replicates was analysed in dupli- 
cate and the average of these two replicates was used for 
further analysis. In case of excessive difference between 
the two replicates (one Ct or more), the run was repeated. 
Considering the high number of samples and genes of 
interest, we opted for this approach in order to reduce the 
number of total runs. Data analysis was performed accord- 
ing to Hellemans et al [46], using the statistical package 
SPSS (IBM). This analysis method takes into account the 
efficiency value of each primer pair. For some genes it was 
necessary to apply a natural log transformation to the 
data, in order to obtain normal distribution of residues. 
To investigate the differences between control and inocu- 
lated samples, we used T-test (p < 0.05). 

Availability of supporting data 

The following files are available on: mynotebook.labarc- 
hives.com. 

Figure 1 - Phylogenetic tree of Rosaceae MLO. 
DOI: 10.6070/H4Z60M0N. 

Additional file 1 - Synteny between apple, peach and 
strawberry. 
DOI: 10.6070/H4TD9V8C. 

Additional files 



experimental replicates. The Ct values have beer^ r^ormalized with three 
reference genes: actin, ubiquitin and elongation factor 1 . Statistical significance 
was determined with a r-test for each individual pair of inoculated and control 
samples at each time point (4, 6, 8 and 24 hpi). The error bars show standard 
errors of the means. Significant differences between inoculated samples and 
control samples are indicated with a *{P < 0.05). 
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